Effect of pairing on one- and two-nucleon transfer below the Coulomb barrier: a 

time-dependent microscopic description. 
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The effect of pairing correlation on transfer reaction below the Coulomb barrier is investigated 
qualitatively and quantitatively using a simplified version of the Time-Dependent Hartree-Fock + 
BCS approach. The effect of particle number symmetry breaking on the description of reaction 
and dedicated methods to extract one and two-nucleon transfer probabilities (Pi„ and P2n) in a 
particle number symmetry breaking approach are discussed. Influence of pairing is systematically 
investigated in the 40 Ca+ 40 > 42 > 44 . 46 > 48 . 50 C a reactions. A strong enhancement of the two-particle 
transfer probabilities due to initial pairing correlations is observed. This enhancement induces an 
increase of the ratio of probabilities P2n/(Pin) 2 compared to the case with no pairing. It is shown 
that this ratio increases strongly as the center of mass energy decreases with a value that could be 
larger than ten in the deep sub-barrier regime. An analysis of the pair transfer sensitivity to the 
type of pairing interaction, namely surface, mixed or volume, used in the theory is made. It is found 
that the pair transfer is globally insensitive to the type of force and mainly depends on the pairing 
interaction strength. 
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I. INTRODUCTION 

The possibility to access cross section much below the 
Coulomb barrier has revealed new aspects like the hin- 
drance of fusion cross section (see for instance [l|) whose 
origin is still debated @, Q. Among possible interpre- 
tation, other competing processes like single- or multi- 
nucleon transfer might eventually be enhanced and/or 
modify the capture process [3, Q. New experimental 
observations [fj-|8[ in the moderate and deep sub-barrier 
regime might lead to important new insight especially on 
the process of pair transfer. The description of such pair 
transfer is particularly complex since it requires to treat 
the quantum tunneling of a composite, eventually corre- 
lated, system. In particular, pairing correlations among 
last bound nucleons is anticipated to play a crucial role. 
Following the pioneering work of Refs [9Hl2|, an impor- 
tant effort is currently being made to improve the de- 
scription of pair transfer in superfluid systems fl3l - l23j ]. 
These approaches have usually in common that transi- 
tion probabilities from the initial to the final nucleus are 
estimated using state of the art Hartree-Fock Bogolyubov 
(HFB) and Quasi-Particle Random Phase Approxima- 
tion (QRPA) nuclear models while the reaction dynam- 
ics part is treated in a completely separated steps using 
coupled channels technique. 

The present work is an attempt to treat nuclear struc- 
ture and nuclear reaction aspects in a common micro- 
scopic framework that includes pairing. Recently, active 
research has been devoted to include pairing correlations 
into the nuclear dynamics using the Time-Dependent 



HFB (TDHFB) approach [15|, |24|, |25j . While current ap- 
plications can be performed in an unrestricted space, due 
to the required effort, applications of TDHFB have been 
essentially made on process involving one nucleus, like gi- 
ant resonances. The use of TDHFB to nuclear reactions 
remains tedious. A simplified version of TDHFB based 
on the BCS approximation is considered. This theory has 
been proposed already some times ago [2?jj and recently 
applied with some success either to collective motion in 
nuclei [23], to reactions in ID models (28J. First step 
toward collisions have been reported in ref. [29| . The 
TDHF+BCS approach has the advantage to be simpler 
than the original TDHFB theory while keeping part of 
the physics of pairing. Note that, time-dependent mi- 
croscopic theories have several advantages compared to 
other techniques. Many effects, like possible dynamical 
deformation or core polarization during the reaction are 
automatically accounted for. In addition, other compet- 
ing phenomena like emission to the continuum and/or fu- 
sion are simultaneously treated. Since, many aspects of 
the theory applied here have been extensively discussed 
in Refs. |27|, [28|], only main aspects features are recalled 
below. 
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II. NUCLEAR REACTIONS WITH PAIRING 



Time-Dependent Hartree-Fock (TDHF) has become a 
standard tool to describe nuclear reactions like fusion 
or transfer reactions (see [3l| and ref. therein). In the 
present work, the TDHF3D code of ref. [HJ is extended 
to include pairing correlations. Below, specific aspects 
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related to the introduction of pairing are discussed. 
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A. Initial conditions 

The reaction is simulated on a 3-dimensional mesh. 
Following the standard procedure [32| , the two nuclei are 
initiated separately and then positioned consistently with 
the desired impact parameter b and center of mass energy 
E c m . The initial wave function can be written as 



|tf(i )) = |$i(*o)>® |$ 2 (*o)> 



(1) 



where |$ a (io)} denotes the many-body wave-function of 
nucleus a = 1,2. Usually, these wave- functions corre- 
sponds to Slater determinants. It is assumed here to 
take the more general form of a quasi-particle vacuum 
written as 

iM*to> = n («fc(*o) + »fc(*o)4(*>)4(*°)) i-)- ( 2 ) 

fc>0 

where o£(to) stands for the creation operator associated 
to the canonical single-particle states, denoted hereafter 
by |<Pfc(£o)) while (ufe(*o)5 u fc(*o)) are the standard upper 
and lower components of the quasi-particle states. Note 
that due to the spatial separation of the two nuclei, a 
common single-particle basis can be used. Accordingly, 
we can omit the a index and directly write the total 
wave-function as: 

i*(*o)> = n ( u fc(*o)+«fc(*o)4(*o)4(*°)) ( s ) 

fc>0 

In practice, initial states for each nucleus have been 
obtained using the EV8 code [33| that solve the self- 
consistent BCS equations in the Energy Density Func- 
tional framework [34| . Single-particle states are written 
in r-space and spin space, denoted by a =7, 4- as: 



[ dr<Pk(r,(r)*U r ), 



(4) 



where \P J. (r) are standard spinors creation operators. 
In EV8, time-reversal symmetry is assumed and single- 
particle states can be grouped by pairs of time-reversed 
states (k, k). Associated quasi-particle creation opera- 
tors are written using the following convention 
for the Bogolyubov transformation: 

A = E r Mr, *o)*j(r) + i>k(r, io)*i(r), 

P\ = Er M*, *o)*l(r) - v k (r, i )*f(r), 
where, using time-reversal properties, we have 

tik (r, i ) = u fc <Pfc(r,t) = Uk<Pk(r,l), 
Ufe(r,t ) = v k <pl(r,l) = u fc ^(r,t). 



(5) 



(6) 
(7) 



The Skyrme Sly4d functional [32| is used in the mean- 
field channel while for pairing, the following effective 
neutron-neutron interaction is used: 

K(r,a;r>') = (l ^^ll±£M^ Sr , r , [1 - P aa , 



where P aa > is the spin exchange operator and where 
Po = 0.16 fm -3 . Here r = n,p stands for neutron or 
proton channel, only neutron-neutron and proton-proton 
interaction are considered. Three different forces, stan- 
dardly called volume (rj = 0), mixed (rj = 0.5) and sur- 
face (ij = 1) will be used below. In each case, the neutron 
pairing interaction strength V™ was adjusted to prop- 
erly reproduce the experimental gap for the calcium iso- 
topic chain deduced from masses using the 5 points for- 
mula [35[ . Theoretical odd systems binding energies have 
been computed using blocking techniques. Values of the 
interaction parameters are reported in ta ble III The pro- 
ton interaction strength is taken from ref. [361 ] but do not 
play any role due to the proton closed shell. Illustration 



interaction 


V 


Vo n [MeV.fm 3 ] 


V PP [MeV.fm 3 ] 


volume 





585 


490 


mixed 


0.5 


798 


755 


surface 


1 


1256 


1462 



TABLE I: Parameters of the neutron-neutron and proton- 
proton pairing strength used in the present work. 



of the pairing gap obtained for the three results of the 
fit are shown in Fig. Q] for the three types of pairing in- 
teraction. A comparison of the neutron pairing gap A« 
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FIG. 1: Experimental (red star) and theoretical neutron gaps 
for Z = 20 as a function of N. The volume, mixed and surface 
interaction are respectively shown by solid, dotted and dashed 
lines. 

obtained using the different interactions is done in Fig. [T] 
The three interactions lead to gaps that are compatible 
with each others and compatible with the experimental 
gaps along the calcium isotopic chains especially in the 
mid-shell. In the following, we consider systematically 
reactions between the doubly magic 40 Ca nucleus and 
other 4x Ca, isotopes. The two reactions 40 Ca+ 40 Ca and 
40 Ca+ 48 Ca will correspond to reactions between two nor- 
mal systems, while in other cases, one of the nucleus will 
present pairing. 
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B. Time-dependent equation of motion 

Once the two nuclei have been initiated, the reaction 
is simulated by performing the dynamical evolution of 
the many-body wave-packet given by Eq. ([3]). Here, the 
TDHF+BCS approximation that may be derived from a 
variational principle [26j or by an approximate reduction 
of the TDHFB equations [23| is used. Since, properties 
as well as numerical aspects of the TDHF+BCS method 
arc discussed in Refs. [13, Hl|, only main ingredients 
of the theory are summarized here. In this theory, the 
wave-function remains at all time in its canonical basis, 
Eq. ([3]) , and the single-particle states evolution identifies 
with the mean-field dynamics with: 



ihd t \(fk) = [h[p] - rjk)\Vk) 



(8) 



where rjk{t) = (tpk(t)\h{p]\Lp k (t)) is a time-dependent 
phase that is conveniently chosen to minimize the effect 
of the £7(1) symmetry breaking. h[p) corresponds here to 
the self-consistent mean-field derived from the Skyrme 
functional including time-odd components. 

Along the dynamical path, the information is con- 
tained in the normal and anomalous densities, denoted 
by p and k written in r-space as: 



'( r ' r ') = ^Z K k{^Pk{r,a)ip- k {r' V) 



fe>0 



-Pfc(»",0>fc(rV))- 



(9) 



(10) 



/ here corresponds to a Fermi distribution with a cutoff 
at 5 MeV and a stiffness parameter equal to 0.5 MeV 
[33|, while 0(rj) equals one for 77 > and zero elsewhere. 
This additional cut-off insures that only states that are 
initially bound are considered during the evolution. 

As discussed in ref. [H|, the reduction of the TDHFB 
to TDHF+BCS leads to some inconsistencies, especially 
regarding the one-body continuity equation, making the 
interpretation of the dynamics difficult. To avoid this 
problem, we used here the Frozen Occupation Approxi- 
mation (FOA). In the FOA, it is assumed that the main 
effect of pairing originates from the initial correlations 
that induce partial occupations of the orbitals and non- 
zero components of the two-body correlation matrix, de- 
noted by C\2- Possible reorganization in time of occu- 
pation numbers and components of C12 arc neglected. 
Said differently, occupation numbers n k and components 
Kfc are kept fixed in time and equal to their initial val- 
ues. Note that similar ideas have been used recently to 
describe two-particle break-up reaction using the Time- 
Dependent Density-Matrix approach (30j. This simpli- 
fication is motivated by the fact that (i) it solves the 
problem of continuity equation [28j (ii) in the simple one 
dimensional model considered in the same reference, it 
gives rather good description of the emission of particles 
and is sometimes more predictive than the full TDHFB 
theory (iii) the FOA approximation applied to collective 
motion in nuclei [13 gives results that are very close to 
the full TDHF+BCS dynamics reported in [27j . 



Illustration of reactions 



(k, k) corresponds to pair of single-particle states that 
were originally degenerated in the static calculation due 
to time-reversal symmetry. n k = V k denote the occu- 
pation numbers while n k = u k v k are the components of 
the anomalous density in the canonical basis. Conjointly 
to the single-particle evolution, the equation of motion of 
the components (u k , v k ) or equiyalently of (n k , K k ) should 
be specified. Following Ref. [27j . we have: 



lh dt nk{t) = ^W A £W-4WA fc (i), (11) 

ihj t n k (t) = K k {t){ m {t) + r^(t)) + A*(t)(2n fc (t) - 1), 

where A k (t) correspond to the pairing field components 
given by 



A fc (t) 



00 



v kkU K k( t )9k(t ) 



(12) 



g k corresponds to the cut-off function that select the pair- 
ing window. This cut-off should be taken consistently 
with the static calculation [33( . Here, a slightly different 
prescription is used compared to the original EV8 with 

9k(t ) = f(Vk(to) - A)/(A - r,k(to))e(- Vk (t )). (13) 



In the present work, we are interested in reactions be- 
low the Fusion barrier like the one presented in Ref. Q 
where the probabilities to transfer x neutrons, denoted 
by P X n can be extracted as a function of the minimal 
distance of approach D during the collision. Assuming 
a Coulomb trajectory, D is related to the center of mass 
energy E'en, through: 



sin(6» cm /2) 



(14) 



where Zp and Zt are the target and projectile proton 
number while 6r m is the center of mass scattering angle. 
Following Ref. [38[ , only central collisions will be consid- 
ered here and different distances D are simulated by vary- 
ing the center of mass energy. Initial conditions are ob- 
tained on a lattice of 2L X x 2L y x 2L Z = 22 A x 22.4 x 22.4 
fm 3 noting that the EV8 code uses symmetries to reduce 
the calculation in one octant of this space. The dynam- 
ical evolution are performed in the center of mass frame 
using a Runge-Kutta 4 algorithm on a spatial grid of 
L x xL y x 2L Z = 60.8 x 22.4 x 22.4 fm 3 with a lattice spac- 
ing Ax = 0.8 fm. The time-step is At = 0.015 x 10~ 22 s. 
Note that, non-equilibrium particle emission is negligible 
due to the small center of mass energy in the entrance 
channel. 
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FIG. 2: Evolution of the neutron density projected onto the 
reaction plan z = 0, for the reaction 46 Ca+ 40 Ca at impact 
parameter b — fm and center of mass energy -E c .m. = 49 
MeV. At initial time (top), t = 20 x l(T 22 s (middle) and 
t — 37 x 10 _22 s (bottom). The neck position is indicated by 
the dashed vertical line. 



As an illustration, the neutron density profiles of the 
reaction 46 Ca+ 40 Ca are shown at different stages of the 
reaction in Fig. [2] During the reaction, the two nuclei 
approach from each other, stick together during a certain 
time and then re-separate. During the contact time that 
strongly depends on the initial center of mass energy, 
they eventually exchange particles. 



D. Particle transfer probability in normal systems 

In practice, the system can be cut into two pieces at 
the neck position to calculate the expectation value of 
the number of exchanged nucleons from one-side to the 
other. By convention, we will denote by B the subspace 
where the lightest nucleus is initially (right side of the 
neck position in Fig. [2]) and B the rest of the total space. 
In a mean-field approach, the simplest way to obtain the 
number of exchanged particles is to estimate the operator 
Nb defined through [39|]: 



(15) 



with the time-dependent wave- function ([3]). Here 0(r) is 
zero on the left side of the neck and 1 elsewhere. 

An illustration of the mean number of transferred par- 
ticle, denoted by N tr = {N B ) - 20, from 46 Ca to 40 Ca 
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FIG. 3: Evolution of the mean number of particles transferred 
from Ca to Ca as a function of time during the reaction 
illustrated in Fig. [2] The probability Pin (dashed line) and 
-p2n (dotted line) to exchange 1 and 2 particles obtained by 
making projection on the side B are also presented as well as 
the quantity P\ n + 2P2 n (open circles). Top: mean- number 
of particles and probabilities obtained without projecting on 
good particle number in the total space. Bottom: same with 
an additional projection on neutron number N — 46 in the 
total space. 



at center of mass energy i? c . m . = 49 MeV is shown in 
Fig. [3] (solid line) . As discussed in ref. HH , a deeper 
understanding of the transfer process can be achieved by 
introducing projection onto good particle numbers in the 
sub-space B (or equivalently B). The projection opera- 
tor on a given number of particles N inside the subspace 
B can be written as (see |38j ) 



Pb{N) 



1 

27 



dtpe 



i<p(N B -N) 



(16) 



where tp is the standard gauge angle. Then the probabil- 
ity to find N particles in the subspace B is: 



P B (N) = (*(t)\P B (N)\*(t)) 

" d V e-^ N (^(t)\^ B (<p,t)) (17) 



1 

JtpN 



where ^s(v?, t)) = e lip B \^>{t)) is a new quasi-particle 
vacuum obtained from the original one by making a ro- 
tation ip in the gauge space from the original state. 

The probabilities extracted by projection arc linked to 
the mean number of particles through the sum rule: 



{N B ) = Y. Np z( N )- 



(18) 



Usually, experimental data are presented in terms of 
probabilities to exchange 1, 2, ... x neutrons (rcsp. pro- 
tons), denoted respectively by P\ ni Pim ■■■ Pxn (resp. 



Pip, P2p, ■■■ Pxp)- In the present reaction, these proba- 
bilities are defined through P xn = -Pb (20 + x) while the 
above sum rule reads N% t = ^2 X xP xn . 

In the present work, probabilities have been evaluated 
using the Pfaffian technique of ref. [4(| and explicit for- 
mulas for the wave-packet are given in appendix El An 
illustration of P\ n and Pi n probabilities obtained using 
the projection method is shown in to p p anel of Fig. [3] for 
the 46 Ca. As already noted in ref. [3g, the In and 2n 
channels are often dominating over other multi-nucleon 
transfer channels leading to iV tr ~ P\ n -\-2Pi n , that is per- 
fectly fulfilled in Fig. [3] after the two nuclei re-separate. 

E. Particle transfer probability in superfluid 

systems 

Strictly speaking, the above method to extract transfer 
probabilities is only valid for normal systems, i.e. when 
the wave-function © identifies with a Slater determi- 
nant that is an eigenstate of particle number. For nuclei 
that present pairing, the initial wave-function explicitly 
breaks the particle number symmetry and the BCS states 
is obtained by imposing the particle number only in aver- 
age. This is for instance the case for the 46 Ca discussed 
above. Said differently, the ground state that is used 
for 46 Ca not only presents a component with N = 26 
neutrons but also with surrounding number of neutrons. 
These components lead to spurious contributions in the 
probabilities extracted in previous section. A possible 
way to remove this contamination is to first select the 
relevant component with No = 20 + 26 particles in the 
full space and then consider the projection onto different 
particle numbers in the sub-space B. In the following, 
we denote by P(Nq) the projector on No particles in the 
full space: 

P(N ) = ±J*d<pe M *- No) , (19) 

where N is now the complete particle number operator. 
More generally, to estimate the possible effect of con- 
tribution from components N =/= No, one can compute 
the probability P(N) that the initial state belongs to the 
Hilbert space of N particles. This probability is defined 
through: 

P(N) = (*(to)|P(JV)|*(t )> (20) 

and is shown in Fig. |4] (top panel). Since 40 Ca has a well 
defined number of particles, by convention, N in the x 
axis of Fig. [4] is taken here as the number of particles of 
its collision partner. Only even components are non-zero 
due to the specific form of the state (Eq. ©). While 
the distribution is properly centered around the imposed 
mean number of particles, non negligible contributions 
coexist, especially for N = N a ± 2 in the initial state. 

To remove possible influence of these spurious compo- 
nents, it is possible to define at all time a state with good 
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FIG. 4: Illustration of the distribution of probabilities to 
have N particles initially for the 40 Ca + 46 Ca reaction in 
the TDHF+BCS case (top) and TDHF with initial filling ap- 
proximation (bottom). Note that since 40 Ca has a good par- 
ticle number, N is defined here as the number of particles in 
the wave-function describing the 46 Ca and is centered around 
N = 26. 

number of neutrons 

\N (t)) ee 1 =P(N )Mt)}. (21) 

y/ {*(t)\P(N )Mt)) 

Then, the mean number of transferred particle as well as 
probabilities P xn can be computed using the same tech- 
nique as in section III Dl Note that, the double projec- 
tion approach proposed here can be regarded as a first 
step towards Projection After Variation (PAV) approach 
standardly applied in nuclear structure generalized here 
to binary reactions. An illustration of result is given 
in bottom panel of Fig. [3] The comparison of projected 
(bottom) and unprojected (top) panel, show that N tr and 
Pin are only slightly affected by the removal of spurious 
components. This is a quite general feature we observed 
in applications presented in the article. However, the dif- 
ference between Pm with and without projection can be 
as large as several orders of magnitude. This conclusion 
also holds for a larger number of particles transferred. 

A second difficulty arises, that could already be seen 
in Fig. [31 when pairing is non-zero. While N^ after colli- 
sions converges to a well defined asymptotic value, small 
oscillations of P\ n and Pi n around their asymptotic val- 
ues remain. These oscillations are also present if the ex- 
pectation value (Ng) is computed as a function of time 
with or without projection onto good particle number in 
the total space. This problem points out a difficulty in 
theories like TDHF+BCS. In a previous article HH, we 
have shown that the one-body continuity equation is al- 
ways respected in TDHFB, while in TDHF+BCS, it is 
respected only if single-particle occupations are frozen, 
which is the case in the present work. However, these 
theories provide only approximate treatment of the two- 
body density matrix and in particular do not respect the 
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two-body continuity equation. This difficulty is not spe- 
cific to the TDHF+BCS theory but is also present in 
TDHFB. Indeed, we have checked in the ID model de- 
veloped in ref. [281 ] . adapted to treat transfer, that sim- 
ilar oscillations occur even if the full TDHFB is solved. 
In the following, results obtained for nuclei with non- 
vanishing pairing will be presented with error-bars with 
height equal to oscillation amplitudes. In most cases dis- 
played below, error-bars will be too small to be seen. 

F. Sensitivity to the pairing residual interaction 

Three different pairing interactions, presented in sec- 
tion III Al have been used to initialize the collision part- 
ners. These interactions lead to different spatial prop- 
erties of the pairing field but have been adjusted to re- 
produce the experimental gaps (see Fig. [T]). In figure 
[5j asymptotic values of one- and two-nucleon transfer 
probabilities are reported as a function of center of mass 
energy for the 40 Ca + 46 Ca for the three pairing inter- 
actions below the Fusion barrier. As seen in the figure, 
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FIG. 5: Comparison of asymptotic probabilities P\ n (solid 
line) and Pi-n (dotted line) obtained with the three pairing 
interactions for the 40 Ca + 46 Ca at various center of mass 
energy below the fusion barrier: volume (open circles), mixed 
(open triangles) and surface (crosses) interactions are repre- 
sented by open circles, open triangles and stars respectively. 
Note that error-bars due to final time oscillations are also 
presented but are very small. 

the extracted transfer probabilities are insensitive to the 
type of interaction used. It turns out, that whatever is 
the form of the pairing effective zero-range vertex, if the 
interaction is carefully adjusted to reproduce the same 
experimental gap (Fig. [IJ, the final transfer rate is also 
the same. Note that the present finding is not in con- 
tradiction with ref. [22[ where different types of inter- 
actions (mixed and surface) were shown to give different 
two-particle transfer from grou nd state to ground state. 
The two forces used in ref. [22[ have been adjusted to re- 
produce the same two neutrons separation energies but 
lead to different pairing gap. In the present work, we do 



not see any evidence of a dependence of the pair transfer 
process on the shape of the pairing force that is used. 

Since all types of force lead to the same probabilities, 
below only results of one of the interaction (mixed) are 
shown. 



G. The no-pairing limit 

Here, we are interested in the enhancement of pair 
transfer probabilities as the pairing is introduced in the 
transport theory. To quantify this enhancement, it is nec- 
essary to also perform calculation without pairing inter- 
action, i.e. TDHF. An additional difficulty arises in the 
comparison between systems with and without pairing. 
Quite often, especially when a given j-shell is partially oc- 
cupied, nuclei initialized with EV8 in the Hartree-Fock 
limit are deformed. The introduction of pairing stabi- 
lizes the spherical shape. Therefore, a direct comparison 
of the case with and without pairing not only probes the 
effect of pairing but also the effect of deformation that is 
(i) not correct for calcium isotopes (ii) not the objective 
of the present work. 

To avoid, possible effects of deformation, we used the 
filling approximation for the last occupied shell, i.e. we 
assume that the last shell has partial occupations rik such 
that all angular momentum projections m are occupied 
in the same way. This insures the convergence of the 
mean-field theory towards non-deformed systems. This 
approach implies that the initial system is not anymore 
described by a wave-packet like in Eq. ((HJ), that would 
identify with a Slater determinant in the usual TDHF, 
but by a many-body density matrix of the form: 

D(t) = lcxp(-^A fe 4(iH(t)) (22) 

fc 

where Z = Tr(cxp(— ^ fc Xka\.(t)ak(t))). The trace here 
is taken on the complete Fock space while a k (t) cor- 
responds to creation operator of the canonical states 
tfk(t). In the filling approximation, the density oper- 
ator corresponds to a statistical density and the infor- 
mation on the system reduces to the knowledge of the 
one-body density matrix p = J2k Wk{t))nk{ipk{t)\ where 
the occupation numbers are related to the coefficients 
through nfc = 1/(1 + e Xk ). The evolution of D(t) is per- 
formed by generalizing the TDHF approach where the 
single-particle states evolve according to the standard 
self-consistent equation of motion (Eq. ([5])) while the 
occupation numbers are kept fixed in time. As far as we 
know, this is the only way to avoid possible mixing of 
deformation and pairing effects and this procedure will 
be taken below as the no-pairing reference. 

Similarly to the pairing case, for non doubly magic nu- 
clei, the density D(t) mixes systems with different parti- 
cle numbers and similar treatment based on double pro- 
jections is necessary to extract transfer probabilities. In 
appendix^ some helpful formulas to perform projection 
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on statistical densities of the form (Eq. (|22|l) are given. 
An illustration of the decomposition of the initial state 
with a mean neutron number (N) = 26 corresponding to 
the 46 Ca is given in bottom panel of Fig. 2J This figure il- 
lustrates that the width of the distribution is comparable 
to the BCS case (top panel) with the difference that odd 
components are also present in the filling approximation. 
Probabilities obtained with the filling approximation will 
be labeled by P xn (MF) while those with pairing will be 
labelled by P xn (BCS). 

As a first illustration of the enhancement of pair trans- 
fer probabilities when pairing is introduced, we have 
extracted systematically the ratios between probabili- 
ties with and without pairing as the pairing interaction 
strength VJj ran is varied in the mixed interaction for the 
reaction 40 Ca+ 46 Ca at E c . m . = 43.7 MeV. These ratios 
are shown in Fig. [6] as a function of V™. When pair- 
ing is accounted for, the two nucleons probabilities have 
been computed using either non-zero components of the 
anomalous density (open triangles) or neglecting them 
(open squares). While the former case corresponds to 
the appropriate treatment of pairing effects, the latter 
case can be regarded as a reference calculation where 
only the sequential transfer of the two neutrons is treated 
while taking properly the occupation number dispersion 
of single-particle states around the Fermi energy. The 
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FIG. 6: Ratio of probabilities P rEn (BCS) / P X „(MF) as a func- 
tion of the pairing strength interaction Vq™™ for the reac- 
tion 40 Ca+ 46 Ca at E c . m . = 43.7 MeV. Pi„(BCS)/Pi n (MF) 
(filled circles) and P 2 „(BCS)/P 2 „(MF) calculated by neglect- 
ing (open squares) or not (open triangles) the anomalous den- 
sity components are shown. 

pairing correlations strongly enhanced the two-particle 
transfer, by an order of magnitude around the physical 
value of the pairing strength (see table H]). Note that the 
enhancement depends on the energy of the collision (see 
below). A smaller but non-zero effect is also seen in the 
one-particle transfer channel. The small increase in P\ n 
stems from the increase of occupation number fragmenta- 
tion as Vq 1 ™ increases. The strong enhancement observed 
when the anomalous density is not neglected compared 
to the case where it is set to zero clearly shows that the 



increase is interpreted as the contribution from direct si- 
multaneous processes. 

III. RESULTS AND DISCUSSION 

In the present work, we have systematically inves- 
tigated the effect of initial pairing correlations on the 
single- and multi-nuclcon transfer by comparing the 
TDHF+BCS with frozen correlations to the mean-field 
dynamics with the filling approximation for collision be- 
tween a 40 Ca and different calcium isotopes below the 
Fusion barrier. In table [TT1 the fusion threshold energy 
B deduced from the mean-field transport theories using 
the technique describes in ref. [4l| are systematically re- 
ported for the different reactions considered here. When 
available, experimental fusion barrier are also shown. It 
is clear from the table, that the introduction of pairing 
has a very weak influence on the barrier height. 



system 


Bo (Exp.) 


Bo (Filling) [MeV] 


Bo (BCS) [MeV] 


40 Ca+ 40 Ca 


53.6 


53.090 


53.090 


40 Ca+ 42 Ca 




52.735 


52.735 


40 Ca+ 44 Ca 


51.8 


52.343 


52.332 


40 Ca+ 46 Ca 




52.069 


52.049 


40 Ca+ 48 Ca 


51.8 


51.935 


51.935 


40 Ca+ 50 Ca 




51.200 


51.247 



TABLE II: Fusion barrier B (in MeV) for the reaction 
40 Ca+ 4:l: Ca. Experimental barrier are taken from the system- 
atic theoretical barrier are computed with a precision of 
0.005 MeV. 



A. Systematic study of two-particle transfer versus 
one-particle transfer 

In FigureEl one- and two-particle transfer probabilities 
obtained for the collision between calcium isotopes are 
displayed as a function of center of mass energy for the 
TDHF+BCS case and no-pairing case. In all cases, when 
one of the collision partner presents pairing, the two- 
particle transfer probabilities are significantly enhanced. 
Conjointly, the one-particle transfer is also increased but 
to a less extend. This implies that the mean number 
of particles exchanged is also influenced by the pair- 
ing correlations due to the sum- rule (Eq. (fT5T) ) . Com- 
paring the TDHF results where the effect of k is in- 
cluded (direct+sequential process) to the case where it 
is neglected (sequential only) , several conclusions can be 
drawn. First, the one-particle probability is almost un- 
changed. Therefore, the enhancement in P\ n observed in 
BCS theory compared to the pure mean-field case is a 
direct consequence of the specific fragmentation of occu- 
pation numbers due to pairing that reduces Pauli block- 
ing effect during the transfer process and is unaffected by 
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36 38 40 42 44 46 48 50 36 38 40 42 44 46 48 50 

E c , m . [MeV] E c . m , [MeV] 

FIG. 7: One (dashed line) and two-particle (solid line) trans- 
fer probabilities as a function of center of mass energy for 
the reactions (a) 40 Ca+ 42 Ca, (b) 40 Ca+ 44 Ca, (c) 40 Ca+ 46 Ca 
and (d) 40 Ca+ 48 Ca. The TDHF+BCS results obtained by ne- 
glecting (crosses) or not (open circles) the anomalous density 
contribution are systematically compared with the mean-field 
case (open triangles). Note that for panel (d), both nuclei are 
closed shell nuclei and pairing correlations vanishes. Accord- 
ingly, only the mean-field result is shown. 

the simultaneous component. A second important con- 
clusion is that the main source of enhancement observed 
in Pi n is coming from the initial two-body correlations 
themselves that lead to direct process during the colli- 
sion. This confirms the observation made in Fig. [BJ 

B. Correlations between two-particle transfer and 
pairing gap 

To further quantify the influence of pairing correla- 
tions on the enhancement of two-particle transfer and 
possible dependence with center of mass energy, the ra- 
tio Pi n (BC S) / Pi n (M F) is displayed as a function of the 
mass of the heaviest nucleus participating to the collision 
and for two different fixed center of mass energies below 
the Coulomb barrier. For comparison, the neutron mean 
gap, 

A Sfc>0 K k&k . , 

&BCS = — , {U) 

obtained for this nucleus is also shown in the top panel. 

Similarly to the pairing gap, this ratio has a typical 
bell shape that drops down to one in magic nuclei. This 
confirms that the enhancement of pair transfer is directly 
proportional to the initial pairing correlations (see for 
instance discussion in ref. [43|). 




40 42 44 46 48 50 

A 



FIG. 8: Top: Mean neutron pairing gap obtained with BCS 
theory for the mixed interaction as a function of mass along 
the isotopic chain. Bottom: Ratio of the two-particle transfer 
probability obtained with and without pairing at fixed center 
of mass energy below the Coulomb barrier reported in table 
HT1 Open triangles and open squares correspond to 4 MeV and 
6 MeV below the Coulomb barrier respectively. 



C. Relationship between Pz n and P?„ 

Experimentally, the no-pairing limit that would be a 
reference for a given reaction, cannot be measured. It is 
therefore important to compare quantities that could be 
measure simultaneously. Usually, the two-particle trans- 
fer Pm is compared to (Pin) 2 , where the latter quantity is 
considered as the probability for a completely sequential 
transfer 0, H, [44|. Such a comparison has the advan- 
tage that both quantities contain all possible effects that 
might influence the transfer of particles as well as pos- 
sible pollution from coming from experimental set-ups. 
In figure [5J this ratio is presented for different theories 
considered here. 

This figure gives interesting insight in the two-particle 
transfer. First, both mean-field and TDHF+BCS where 
only the fragmentation of single-particle state is ac- 
counted for while Cn = 0, lead to almost identical ratios. 
This aspect was not clear from Fig. [7] where different 
fragmentations obtained with the filling approximation 
and from BCS with C\i = lead to differences for both 
Pin and Pin- The mean- field theory or equivalcntly the 
BCS where initial correlations are neglected could be con- 
sidered as a way to mimic independent transfer of the 
two-particles. 

It turns out that simple combinatorial arguments can 
be used to understand analytically the sequential limit. 
Let us denote by p the average probability to transfer one 
particle from the 4a: Ca to 40 Ca. Here "average" means 
that we disregard the fact the the probability depends 
on the initial and final single-particle states. It turns out 
that the total probability to transfer 1, 2, k nuclcons 
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FIG. 9: Ratio P 2 n/(Pm) 2 as a function of beam ener- 
gies. The different panel correspond to different reactions: 
(a) 40 Ca+ 42 Ca, (b) 40 Ca+ 44 Ca, (c) 40 Ca+ 46 Ca and (d) 
40 Ca+ 48 Ca. The TDHF+BCS results obtained by neglecting 
(crosses) or not (open circles) the anomalous density contri- 
bution are systematically compared with the mean-field case 
(open triangles). 



during the reaction Ca+ Ca for x > 2 in the MF 
approximation is consistent with: 



Pin = ^inP(l - p) 

Pin = ^2nP 2 (l ~ p) 



N v -1 



AL-2 



n knP k (i- P ) 



7V„-fe 



(24) 



where N v = x is the number of valence nucleons (with the 
constraint k < N v ) in the emitter with respect to the in- 
ert core of Ca, while fifc n is a purely combinatorial fac- 
tor that depends on the number of nucleons in the valence 
shell and on the number of available single-particle states 
in the / 7 / 2 empty shell of the receiver nucleus (Nf = 8). 
f2fc„ simply counts the number of possibilities to select 
k particles among N v times the number of ways to put 
them in the /' 1 1 2 shell, i.e. 



Ay 



A; ■/■) 



xN f (N f -l)---(N f -k + l). 



k\(N v - k)\ 
Accordingly, one can anticipate that 

P2n 1 (N v - 1) (Nf - 1) 

(P ln f 2 N v N f " (1-p) 

^ 1 (N v - 1) (Nf - 1) 
2 N v Nf 

where the last approximation holds if p <C 1. 



(25) 



This simple approximation turns out to work very well 
in the mean-field case (or equivalently in the pairing 
case when k is neglected). In figure [TUl the quantity 
P2n/(Pin) 2 is compared to the left side of Eq. (f2~5j) for 
the different reactions considered here. We see that for 
a wide range of center of mass energy, mean-field results 
perfectly matches the relation (j2"5j) . The fact that such a 
simple description is adequate in mean-field theory is not 
trivial. Indeed, in this theory, nucleons are quantal ob- 
jects interacting first with two cores (the emitter and the 
receiver nucleus) that are not fully inert and second with 
each other through the self-consistent mean-field. Last, 
the two transferred nuclei are fermions and are subject to 
the Pauli exclusion principle. This induces automatically 
correlations during the transfer. If a particle is already 
transferred to a certain single-particle level, this auto- 
matically forbid the other particles to be transferred to 
the same level. The latter effect is automatically included 
in the present theory and partially described through the 
factor fi fen in Eq. 




42 44 46 
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10: Ratio P 2 „/(Pm) 2 

42 Ca+ 40 Ca (open circles) 44 



as a function of center of mass 
Ca+ 40 Ca 



FIG. ru: naxio r 2n 
energy for the reaction 
(open triangles), 46 Ca+ 40 Ca (cross) 
squares). The horizontal lines correspond in each case to the 
value of left side of Eq. (|25[) where Nf = 8 while N v = x for 
4x Ca-|- 40 Ca reactions. 



and 4S Ca+ 4U Ca (open 



Focussing now on the results accounting for initial cor- 
relations (open circles in Fig. [S]), a strong, center of mass 
energy dependent, enhancement of the ratio is seen. The 
ratio increases significantly as the energy decreases from 
a value lower than 1 up to 20 in some cases. The present 
enhancement is at variance with the recent experimental 
observation in 40 Ca+ 96 Zr where Pi n ~ 3(Pi„) 2 has been 
observed almost independently of the center of mass en- 
ergy 0- It is worth mentioning however that the one- 
and two-particle transfer is anticipated to depend signifi- 
cantly on the structure properties, single-particle energies 
and spectroscopic factors, of the two collisions partners. 

In addition, here we are focussing on pairing correla- 
tion effect and paid a particular attention to not mix 
effects coming from static deformation in nuclei. Last, 
mean-field alone cannot grasp the physics of the quan- 
tum fluctuations in collective space. The inclusion of 
pairing partially cure this problem by increasing fluctu- 
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ations of two-body observablcs. However, pairing alone 
does not contain all physical effects to treat this problem. 
This is clearly illustrated close to magicity where pairing 
vanishes. In that case, TDHF dynamics is known to fail 
to reproduce transfer cross section. Recently, a stochas- 
tic mean-field approach was shown to properly describe 
quantal collective fluctuations especially in magic nuclei 
[45l - |48j and leads to realistic description of the nucleon 
exchange process. It would be interesting, in the near 
future to explore the possibility to combine stochastic 
methods with the present BCS approach. 

IV. CONCLUSION 

The TDHF+BCS theory with frozen correlations is 
used here to investigate the effect of pairing on onc- 
and two-nucleon transfer below the Coulomb barrier. A 
method based on projection onto particle number is de- 
veloped to properly extract transfer probabilities from 
theories that break the U(l) symmetry. In addition, a 
particular attention is paid to compare with a no-pairing 
limit free from possible effect of deformation. With this 
technique, the enhancement of two-particle transfer due 
to pairing correlations is studied qualitatively and quanti- 
tatively for reactions involving different calcium isotopes. 
It is shown, that when one of the collision partner has 
non-zero pairing, a strong enhancement of pair transfer 
is observed. This increase is directly proportional to the 
initial pairing correlations in the superfiuid nucleus and 
turns out to strongly depend on the center of mass en- 
ergy 

Appendix A: Formulas for projection 

In the present appendix, formulas useful for the nu- 
merical estimate of particle number projection are given 
for many-body quasi-particle states and density opera- 
tors respectively given by Eqs. ([3]) and (f2"2"j). 

1. Particle number projection of density operators 

Starting from the density (|22p , the probability to have 
N particles in the subspace B can be written as: 

Pb ^ = "hj di P e ' itpNrTr { eiipflB D) ■ ( ai ) 

The estimate of the trace can be made by writing the 
operator Nb in the canonical basis {fi} associated to 
the density. Using the expression of Nb and the fact 
that the canonical basis forms a complete basis of the 
total single-particle space, it could be easily shown that: 

(A2) 



where 

°fi = Hl ^*(r,<7)^(r,a)6(r) = (i\j) B . (A3) 

Then, using formula (A. 16) of ref. [49| leads to 

Tr(e^ B £») = - exp[Trln(l + e** 3 * e~ M )] 

= -det(l + e-^° B e~ M ) (A4) 
z 

where 

(e-^° B ) = Faifp) = % + (i\j) B (e^ - 1), (A5) 
while from formula (8.11) of ref. [49|, we have: 

(e~ M h = ^~ (A6) 

1 Tli 

and z = n^(l + T^r - )• Altogether, we obtain: 

Pb(N) = A- [*dw-* N tet(0--n j )6 ii +F ij (<p)n j ). 
^ Jo 

Note that in the case where the statistical density iden- 
tifies with a Slater determinant (n, = 0, 1), the formula 
given in ref. [38j is properly recovered. Formulas for 
the double-projection technique can be derived using a 
similar technique. 

2. Projection with quasi-particle states 

To perform projection of quasi-particle vacuum onto 
good particle number, we used the recently proposed 
Pfaffian method [4(], Ho, [H| • Since the Pfaffian technique 
has been largely discussed recently, here, only specific for- 
mulas useful in the present article are given. Again, we 
first consider the projection on the B subspace as an il- 
lustration. We need to perform the overlap between the 
quasi-particle state ([3]) and its gauge angle rotated coun- 
terpart: 

i*) = n («*+«*44) 

fe>0 

l**(¥0) = l[(u k + v k b{(v)bl(<p))\-), 

k>0 

where 

b \^) = E / dve^ip^a^Hx), (A7) 
= £>^)at. (A8) 

3 

The matrix F plays the role of the matrix R in ref. 
and the overlap between the non-rotated and rotated 
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state are given by Eq. (5) of this reference. In the present 
case, we obtain: 



271 



4^' ^ 



2n 



dif'e 



-iipN-<p'N' 



(*o|*bM> 



-l) n 



n 



y2 



pf 



K. M(ip) 
-M\<p) -K,* 



where /C and M. are matrix of size 2n x In where n is 
the number of single-particle states with i > 0. These 
matrices can be decomposed in 2 x 2 matrix blocks as: 



K. = 







and 



M{<p) 



[vp>jFij(.<p)] bnVjFijiv)] 



where matrix elements are directly indicated in each nxn 
block. 

For the double projection, the probability to find N' 
particles in the space B for a system of N particles in 
the total space is given by 

P b{ N,N>) = (iV| ^ (7V/)|7V) = mPB(N')P(N)\^ } 



(N\N) 



(V\P{N)\V) 



Therefore, we need to integrate with respect to two gauge 
angles. 

(■$\P B (N')P(N)\iZ>) = 



where (^f\^>B(f,f')) can be calculated using formula 
(|A9|) except that Fij(<p) is now replaced by Fij(<p,(p') = 

Numerically, the gauge integral are discretized using 
the Fomenko method [52| with 20 points. Note that 
during the time evolution, due to accumulated numeri- 
cal errors, small violation of orthonormalization between 
single-particle states can occur, this might lead to large 
errors in the extracted transfer probabilities. To avoid 
this problem, a Gram-Schmidt orthonormalization algo- 
rithm is used prior to applying the Pfaffian formula. 
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